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LONG-TERM  GOALS 

The  goal  of  our  Tidal  Flats  research  is  to  expand  our  understanding  of  the  erosional  and  depositional 
processes  that  lead  to  exchange  of  mud  between  tidal  channels  and  tidal  flats  and  its  impact  on 
sediment  strength. 

OBJECTIVES 

Our  research  on  Tidal  Flats  DRI  has  two  primary  objectives: 

1 .  Improve  understanding  of  how  turbulence  and  sediment  concentration  affect  the  size-dependent 
depositional  flux  of  sediment  to  the  seabed  on  tidal  flats. 

2.  Improve  understanding  of  how  sediment  texture  and  sorting  in  the  seabed  affect  the  size-dependent 
erosional  flux  from  the  seabed  on  tidal  flats. 

APPROACH 

We  have  focused  our  studies  on  two  compositionally  different  sites  at  the  Southern  end  of  Willapa 
Bay,  WA,  that  were  identified  during  a  regional  grain  size  survey  conducted  in  September  2008 
(Figure  1).  The  first  is  a  muddy  site  located  near  the  secondary  channels  and  flats  west  of  the  Bear 
River  channel  near  Round  Island.  The  second,  sandier,  more  exposed  site  is  located  at  a  lower  tidal 
level  approximately  2.5km  to  the  northwest  of  the  first  site.  Cores  for  erodibility  studies  were 
collected  from  the  channels  and  flats  from  both  sites.  Suspended  sediment  measurements  have  been 
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concentrated  at  the  muddy  site  in  order  to  collect  spring-neap  time  series  and  to  help  understand 
changes  in  flux  that  may  result  in  low  strength  muds. 


Figure  1:  Location  map  for  Willapa  Bay  Tidal  Flats  coring  and  suspended  sediment  observations. 

Yellow  denotes  areas  where  cores  were  collected  for  erodibility  studies  in  September  2008  and 
March  and  July  2009.  Red  shows  area  where  suspended  sediment  measurements  and  erodibility 
studies  were  carried  out  in  March  and  July  2009.  Labels  indicate  geometric  mean  diameter  of 
disaggregated  inorganic  bottom  sediments  and  volume  percentage  less  than  4  pm  in  diameter.  The 
site  SSW  of  Round  Island  is  finer  grained  than  the  site  NW  of  Round  Island.  Image  is  ~4km  wide. 

WORK  COMPLETED 

To  examine  seasonal  variations  in  sediment  dynamics,  experiments  were  carried  out  on  the  Willapa 
tidal  flat  in  September  2008  and  March  and  July  2009.  Sediment  cores  were  collected  to  determine 
erodibility  at  stations  located  on  both  the  tidal  flat  and  in  channels.  Cores  with  undisturbed  sediment- 
water  interfaces  were  obtained  using  a  specially  designed  hand  held  corer  from  small  boats  and  kayaks 
A  Gust  erosion  chamber  was  used  to  erode  sediment  at  known  shear  stresses  and  this  material  was 
collected  for  disaggregated  inorganic  grain  size  (DIGS)  analysis.  In  addition  a  regional  survey  of 
surficial  sediment  grain  size  was  conducted  in  September  and  July.  In  September  and  March  a  LISST 
lOOx,  digital  floe  camera  (DFC),  and  a  digital  video  camera  (DVC)  attached  to  a  settling  column  were 
deployed  to  examine  in-situ  particle  size  and  settling  velocity  (Figure  2).  Measurements  were  made  at 
both  the  sand  and  mud  flat  sites.  In  July,  a  second  LISST  and  DFC  were  added  to  the  survey  to  allow 
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simultaneous  measurements  of  in-situ  particle  size  in  both  the  channel  and  flat  at  the  mud  site  and  to 
make  simultaneous  measurements  on  the  flat  at  the  sand  and  mud  sites. 


Figure  2.  Instruments  on  the  Willapa  Bay  tidal  flat.  From  left  to  right,  Digital  Settling  Velocity 

TM 

Camera  (DVC),  LISST  lOOx,  Digital  Floe  Camera  (DFC),  and  2  Nortek  Aquadopps  (on  loan  from 

Washington  State  University  in  Vancouver). 


RESULTS 

Several  basic  results  are  emerging  from  our  work  to  date. 

•  Suspended  sediment  concentration,  floe  size,  and  settling  velocity  are  maximal  when  currents 
are  strongest.  This  observation  indicates  that  gross  downward  flux  due  to  sinking  is  maximal  at 
peak  currents  (Figure  3). 

•  Suspended  sediment  concentration,  floe  size,  and  settling  velocity  are  larger  on  flood  tides  than 
on  ebb  tides.  This  observation  indicates  that  there  is  asymmetry  in  vertical  and  horizontal  flux 
at  the  height  of  our  sensors  (0.5  mab),  with  net  flux  of  fine  sediment  onto  the  tidal  flats  (Figure 
3). 

•  Size  sorting  during  erosion  depends  on  bed  sediment  texture.  When  the  percentage  of  sediment 
smaller  than  4  pm  is  less  than  5-10%,  sediments  behave  non-cohesively.  When  this  fraction  is 
higher,  all  sediment  smaller  than  ~16  pm  is  resuspended  equally  (Figure  4;  Law  et  ah,  2008). 
This  observation  indicates  that  sandy  sites  will  remain  sandy,  and  muddy  sites  will  accumulate 
more  mud. 
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•  The  DIGS  for  surficial  sediment  samples  collected  in  March  has  a  larger  proportion  of  floe 
settled  material  than  samples  from  September.  July  samples  still  require  analysis  to  determine 
trend  on  a  longer  temporal  scale. 
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Figure  3.  Time  series  of  water  level,  beam  attenuation  and  particle  groups  on  a  muddy  tidal  flat  in 
Willapa  Bay.  In  the  top  panel,  as  water  level  (black  line)  increases,  beam  attenuation(gray  line) 
increases.  At  high  water  beam  attenuation  decreases  indicating  that  sediment  is  settling  from 
suspension.  As  water  level  falls  on  the  ebb,  attenuation  increases  again,  but  not  to  the  magnitude 
that  it  had  on  the  flood  tide.  In  the  second  panel,  a  time  series  of  bottom  current  speed  is  plotted 
using  a  fill  color  that  corresponds  to  the  dominant  particle  size  entropy  group  shown  in  the  bottom 
panels.  Entropy  analysis  groups  particle  size  spectra  in  a  way  that  maximizes  both  similarity  among 
size  distributions  within  a  group  and  differences  among  size  distributions  between  groups 
(Mikkelsen  et  al.,  2007).  Particle  size  distributions  indicate  that  floes  are  larger  on  the  flood  than  on 
the  ebb.  The  asymmetry  between  beam  attenuation  and  particle  size  on  the  flood  versus  the  ebb  can 
not  be  explained  by  shoreward  pumping  of fine  sediment.  It  may  indicate  that  there  is  flow 
stratification  on  the  ebb  tide  that  prevents  floes  from  diffusing  to  the  height  of  the  particle  sensors 

0.5  mab. 

The  observations  fit  into  a  conventional  model  of  sedimentation  on  mudflats.  According  to  the 
conventional  model,  sands  on  tidal  flats  are  found  below  the  mid-water  level,  on  parts  of  the  seabed 
that  experience  maximum  tidal  flows.  The  sand  site  and  the  channels  at  the  mud  site  are  in  this 
category  (Figure  1).  Muds  are  found  higher  on  the  flats.  Because  of  the  longer  period  required  for 
muds  to  sink  from  suspension,  they  are  carried  landward  on  flooding  tides,  and  they  deposit  at  high 
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slack  water.  Because  of  cohesion,  freshly  deposited  muds  are  not  re-eroded  on  the  ebb  tide.  This 
model  generally  fits  our  observations  of  particle  size  and  size  sorting  during  erosion.  The  hypothesis 
that  freshly  deposited  mud  on  the  upper  flats  is  not  re-eroded  on  ebb  tides,  however,  is  difficult  to 
reconcile  with  other  observations. 

In  general,  sediment  on  the  muddy  flats  is  difficult  to  erode,  suggesting  that  there  is  little  freshly 
deposited  sediment  at  the  surface.  In  contrast,  sediment  in  the  channels  tends  to  be  more  mobile.  This 
observation  could  be  due  to  the  cohesive  behavior  of  sediment  on  the  muddy  flats  versus  the  non- 
cohesive  behavior  at  sandier  sites.  This  hypothesis  requires  rapid  and  significant  cohesion  of  freshly 
deposited  muds,  and  it  has  difficulty  explaining  the  presence  of  easily  eroded  muds  in  channels  and  at 
sandy  sites.  According  to  the  conventional  model,  the  easily  eroded  muds  should  be  exported  to 
muddy  flats  and  remain  there,  which  in  the  absence  of  a  significant  sediment  source  would  leave  the 
sandy  flats  and  channels  depleted  of  fine  mobile  sediment.  An  alternative  hypothesis,  supported  by 
Chuck  Nittrouer’s  observations  from  a  jack-up  barge,  is  that  freshly  deposited  mud  on  the  upper  flats  is 
re-eroded  on  ebb  tides  and  is  returned  to  the  lower  flats  and  channels.  This  suspended  material  is  not 
measured  by  our  particle  sensors  0.5  mab  because  of  water  column  stratification,  which  limits  the 
upward  turbulent  diffusion  of  floes  and/or  generates  differential  floe  breakage  on  flood  versus  ebb 
tides  because  of  suppression  of  turbulence  (e.g.,  Scully  and  Friedrichs,  2007). 

Resolving  the  fate  of  fine  sediment  advected  onto  the  upper  flats  is  critical  to  the  prediction  of 
trafficability.  If  freshly  deposited  fine  sediment  remains  on  the  upper  flats,  then  these  areas  pose  the 
greatest  problem  for  trafficability.  If,  on  the  other  hand,  fine  sediments  return  to  the  channels  on  each 
ebb  tide,  then  these  areas  are  more  prone  to  accumulation  of  low-trafficability  muds.  We  hypothesize 
the  latter  scenario,  and  future  work  will  test  this  hypothesis  by  conducting  simultaneous  measurements 
of  sediment  flux  and  erodibility  on  the  flats  and  in  an  adjoining  channel. 

IMPACT/APPLICATION 

Physical  properties  of  the  seabed  depend  fundamentally  on  particle  size.  This  dependence  endows 
energetic  coastal  tidal  flats  with  extraordinary  variability  in  seabed  properties,  because  these  are  the 
only  environments  on  Earth  where  mud,  sand,  and  gravel  are  all  actively  eroded,  transported,  and 
deposited.  Over  the  course  of  a  single  tide,  a  spring-neap  tidal  cycle,  a  season,  or  a  year,  seabed  grain 
size  at  a  single  location  can  change  dramatically  on  an  energetic  flat  (e.g.,  Choi  et  ah,  2004). 

Variability  in  seabed  properties  represents  a  major  concern  for  safe  and  efficient  travel  across  tidal 
flats.  Especially  problematic  are  “oozing  mudflats”  (cf.  Wells,  1983,  Frey  et  ah,  1989)  with  high 
porosity  and  low  strength.  People  and  equipment  can  sink  deeply  into  such  muds,  making  travel 
difficult  and  potentially  dangerous.  The  link  between  seabed  properties  and  “trafficability”  motivates 
efforts  to  understand  the  formation  of  low-strength  muds  so  that  models  that  can  predict  the  locations 
of  such  deposits  can  be  developed. 
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Figure  2.  Results  from  erosional  sorting  studies  carried  out  on  three  cores  from  Willapa  Bay  in 
September  2008  collected  at  three  sites:  top  =  sand,  middle  =  mud  flat,  bottom  =  mud  channel.  Left 
panels  show  DIGS  of  bottom  sediment  (black)  and  material  eroded  at  increasing  stress  (colored) 
with  corresponding  geometric  mean  diameter  percentage  of  sediment  <  4  pm.  Right  panels  show, 
size-specific  relative  mobilities  (SSRM)  for  the  eroded  cores.  The  SSRMs  divide  the  relative 
concentration  in  a  size  class  in  suspension  by  its  relative  concentration  in  the  bed.  Values  greater 
than  one  indicate  preferential  erosion  of  that  size  class,  while  a  value  of  one  indicates  that  the 
particle  size  has  the  same  relative  concentration  in  suspension  as  it  does  in  the  bed,  i.e.  there  has 
been  no  size  sorting.  The  pattern  in  the  sand  core  is  consistent  with  non-cohesive  behavior,  which 
generally  occurs  when  the  fine  fraction  is  less  than  5-10%  (Law  et  al.,  2008).  The  channel  core 
shows  SSRMs  near  to  unity  over  a  range  of  size  classes,  indicating  cohesive  behavior  and  minimal 
size  sorting.  This  behavior  is  consistent  with  its  large  fine  fraction.  The  GUST  1  core  shows 
intermediate  sorting  behavior,  consistent  with  its  fine  fraction  of  9%. 
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RELATED  WORK 


This  work  is  being  carried  out  in  close  collaboration  with  other  ONR  sponsored  researchers.  Cores 
collected  were  also  measured  for  erosion  rates  and  cumulative  mass  eroded  (Pat  Wiberg,  UVA), 
porosity  and  organic  content  (Rob  Wheatcroft,  OSU),  and  fracture  strength  (Bernie  Boudreau  and 
Bruce  Johnson,  Dalhousie).  Chuck  Nittrouer  carried  out  a  regional  coring  program  to  determine 
sediment  source  and  accumulation  rates  on  the  Willapa  flats.  Results  from  the  in-situ  measurement  of 
particle  size,  transmission  (Cp),  and  settling  velocity  are  being  applied  to  the  ONR  funded  Ocean 
Optics  OASIS  project. 
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